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We study the effect of impurities on the ground state and the iow-temperature Ohmic dc transport 
in a one-dimensionai chain and quasi-one-dimensional systems of many parallel chains. We assume 
that strong interactions impose a short-range periodicicity of the electron positions. The long- 
range order of such an electron crystal (or equivalently, a 4fc_F charge-density wave) is destroyed 
by impurities, which act as strong pinning centers. We show that a three-dimensional array of 
chains behaves differently at large and at small impurity concentrations N. At large N, impurities 
divide the chains into metallic rods. Additions or removal of electrons from such rods correspond to 
charge excitations whose density of states exhibits a quadratic Coulomb gap. At low temperatures 
the conductivity is due to the variable-range hopping of electrons between the rods. It obeys the 
Efros- Shklovskii (ES) law — ln<r ~ (Tes/T) 1 '' 2 . Tes decreases as N decreases, which leads to an 
exponential growth of a. When N is small, the metallic-rod (also known as "interrupted-strand") 
picture of the ground state survives only in the form of rare clusters of atypically short rods. They 
are the source of low-energy charge excitations. In the bulk of the crystal the charge excitations 
are gapped and the electron crystal is pinned collectively. A strongly anisotropic screening of the 
Coulomb potential produces an unconventional linear in energy Coulomb gap and a new law of the 
variable-range hopping conductivity — In it ~ (Ti/T) 2 ^ 5 . The parameter Ti remains constant over 
a finite range of impurity concentrations. At smaller N the 2/5-law is replaced by the Mott law, 
-In cr ~ (T M /T) 1/4 . In the Mott regime the conductivity gets suppressed as N goes down. Thus, 
the overall dependence of a on N is nonmonotonic. In the case of a single chain, the metallic-rod 
picture applies at all N. The low-temperature conductivity obeys the ES law, with log-corrections, 
and decreases exponentially with N. Our theory provides a qualitative explanation for the transport 
properties of organic charge-density wave compounds of TCNQ family. 



I. INTRODUCTION 

In recent years electron transport in quasi-one- 
dimensional (quasi-lD) systems moved into focus of both 
fundamental and applied research. Quantum wires, nan- 
otube ropes, conducting molecules, etc., are being ex- 
amined as possible elements of miniature electronics 
devices. In parallel, discovery of quasi-lD structures 
termed "stripes" in correlated electron systems (high-T c 
cuprates, quantum Hall devices, etc), invigorates efforts 
to understand unconventional phases in two and three 
dimensions starting from models of weakly coupled ID 
chains. 

Experimentally, the low-temperature conductivity 
cr(T) of quasi-lD systems is often of the insulating type. 
Its temperature dependence gives information about the 
nature of charge excitations. For example, the activated 
dependence — ln<j(T) oc 1/T indicates a gap in the spec- 
trum. In quasi-lD systems such a gap commonly arises 
from the Mott-Peierls mechanismjAiii where the commen- 
surability with the host lattice is crucial. Yet there are 
many ID and quasi-lD systems where commensurability 
plays a negligible role. In this situation the jellium model 
(an electron gas on a positive compensating background) 
is a good approximation. This is the kind of systems 
we study in this paper. We will show that their low- 
temperature transport is dominated by a variable-range 
hopping (VRH), which leads to a slower than exponen- 



tial T-dependence2i4 of the conductivity. Our theory may 
apply to a number of systems, both naturally occuring 
and man-made. Prototypical ID examples are individual 
quantum wires an carbon nanotubes. Stripe phases^*& 
quantum wire arrays in heterojunctions^ carbon nan- 
otube films' and atomic wires on silicon surface^ are 
two-dimensional (2D) examples. In three dimensions 
(3D), an important and well studied class of quasi-lD 
compounds are charge-density wavesji2*ii*i£ (CDW). 

To characterize the strength of Coulomb correlations in 
a quasi-lD system we define the dimensionless parameter 
r s = a/2as, where a is the average distance between elec- 
trons along the chain direction and as = ft 2 n/me 2 is the 
effective Bohr radius. The latter is expressed in terms of 
the dielectric constant of the medium n and the electron 
band mass m. In practically all known realizations of ID 
and quasi-lD systems, r s exceeds unity, often by orders 
of magnitude. Below we assume that r s 1. Under this 
condition the dynamics of electrons can be treated semi- 
classically. Neglecting quantum fluctuations altogether 
for a moment, we arrive at the picture of electrons form- 
ing a classical ID Wigner crystal in the case of a single 
chain (Fig. 2Jl) or an array of such crystals in 2D and 3D 
(Fig. [iy> and c) , with the period along the chain equal 
to a. Formation of the crystal enables the electrons to 
minimize the energy of their mutual Coulomb repulsion. 
In order to correctly assess the role of quantum fluctua- 
tions (zero-point motion) in such crystals, one has to take 



2 



into account two circumstances. First is the unavoidable 
presence of random impurities that act as pinning cen- 
ters. Second is a finite interchain coupling (in 2D and 3D 
systems). Either one is sufficient to make the quantum 
fluctuations of electron positions bounded. For exam- 
ple, in the case of a single chain, the slow growth of the 
zero-point motion amplitude with distanc e i 1 ^ is termi- 
nated at the nearest strong pinning center. In higher di- 
mensions, the zero-point motion amplitude is finite even 
without impurities because of the interchain interaction. 
In fact, renormalization group approaches^ indicate that 
the Wigner crystal or, equivalently, 4fci?-CDW is the true 
ground state of a system of weakly coupled chains start- 
ing already from rather modest r s . In all situations, the 
net effect of quantum fluctuations is to slightly renor- 
malize the the bare impurity strength and/or bare in- 
terchain coupling. The calculation of renormalized pa- 
rameters is possible via the standard bosonization tech- 
nique^ For the treatment of, e.g., ID case, one can con- 
sult Refs. Below such a renormalization is 
assumed to be taken into account and it is not discussed 
further. Henceforth we will often refer to the systems we 
study as electron crystals. 

Due to impurity pinning, at zero temperature and in 
the limit of small electric field (Ohmic regime) the con- 
ductivity of the electron crystal vanishes. This behavior 
is common for all pinned systems. It motivated a large 
body of workiSii^ devoted to mechanisms of nonlinear 
transport that become possible in the presence of a fi- 
nite electric field, in particular, creep and sliding. In the 
present context, such mechanisms would involve a collec- 
tive motion of large numbers of electrons that overcome 
pinning barriers either by thermal activation or by quan- 
tum tunneling. It has been understood, however, on some 
qualitative level^ that if compact charge excitations are 
allowed by the topology of the system, then such excita- 
tions would dominate the response at low temperatures 
and would give rise to a nonzero Ohmic conductivity at 
T > 0. Below we will demonstrate that this is indeed 
the case in the electron crystals. We clarify the nature of 
the compact low-energy charge excitations and propose 
a theory of their low-temperature Ohmic transport that 
consistently addresses the role of long-range Coulomb in- 
teractions |Si 

In this paper we are focused exclusively on the charge 
transport and ignore any effects related to the spin degree 
of freedom. This is legitimate for r s 3> 1 because electron 
are tightly localized at the sites of the classical Wigner 
crystal and the energy of their spin-dependent exchange 
interaction is exponentially small. 

We will assume that impurities that pin the crystal 
are strong enough to enforce preferred order of electrons 
nearby or, in the CDW terminology, the preferred phase. 
The relation between the phase cf> and the elastic displace- 
ment of the crystal u is <f) = — {2n/a)u. In Fig.2]impu- 
rities are shown by vertical tick marks and it is assumed 
that they interact with nearby electrons by a repulsive 
potential comparable in magnitude to the Coulomb in- 
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FIG. 1: Pinned ID and quasi-lD systems on a uniformly 
charged positive background. Dots and tick marks label the 
positions of electrons and impurities, respectively, (a) ID 
crystal. The + and — signs denote the charge of the metal- 
lic rods, which have an average length of I. (b) An array of 
decoupled chains in the case I < l 3 - (c) An array of cou- 
pled chains for the case I > l s - The preferred arrangement 
of the electrons on neigborng chains may depend on the ex- 
act geometry of the system. In order not to complicate the 
drawing, we adopted the convention where the ground state 
corresponds to the same horizontal positions of electrons on 
each chain. The interchain interactions try to diminish the de- 
viations from this ground state leading to dipolar distortions 
of a characteristic length l 3 around impurities. 



teraction energy e 2 /na between nearest electrons on the 
chain. This condition is sufficient to ensure that the im- 
purity to act as a strong pinning center. An example of 
such an impurity is an acceptor residing on the chain. In 
the ground state one electron is bound to the acceptor 
and the electron-acceptor complex (of total charge e) is 
built into the crystal, i.e., it is positioned squarely in be- 
tween the two closest other electrons. One can say that 
the crystal contains a plastic deformation — a vacancy 
bound to the negatively-charged acceptor. Overall, the 
region around the impurity is electrically neutral. 

In the case of a single chain (Fig. ^l), strong impu- 
rities divide the crystal into segments, which behave as 
individual metallic rods. A charge can easily spread over 
the length of each rod, while it has to tunnel through an 
impurity to move to a neighboring rod. Each rod con- 
tains an integer number of electrons but the charge of the 
positive background is random because of the assumed 
incommensurability. In the ground state the distribution 
of the rods' total charges, electrons plus background, is 
uniform between — e/2 and +e/2 (</> is between —tt and 
7r). Larger charges cost more Coulomb energy and cor- 
respond to charge excitations above the ground state. 
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Transitions between ground and excited states occur by 
discrete changes in the number of electrons on the rods. 

Let e be the minimal by absolute value change in the 
self-energy of a given rod due to change of its charge 
by one unit. By the self-energy we mean the Coulomb 
energy of interaction among the electrons on the given 
rod, charges on all other rods held fixed. We adopt the 
convention that the chemical potential (Fermi energy) 
corresponds to the zero energy. In this case, e is non- 
negative (nonpositive) for addition (subtraction) of the 
electron. We denote by 5s (e) the distribution function of 
e averaged over impurity positions. We refer to gs as the 
bare density of states of charge excitations. In Sec. IIIII 
we show that random distribution of charges of the rods 
creates a finite gs (0). Small e come from rods with net 
charges close to ±e/2. These rods make possible a VRH 
conductivity at low temperatures. 

Consider now a 3D system of parallel chains. Impuri- 
ties with concentration N divide the chains into segments 
of average length I = 1/Na\, where a\ is the area per 
chain in y-z plane (we assume the chains to be along the 
x-direction) . We get two cases distinguished by the rel- 
ative importance of interchain interactions. In the first 
case (Fig.njj), the chains are far enough from each other 
and/or the impurity concentration is large enough so that 
the interchain coupling over the length ~ I of a typical 
segment can be neglected in comparison with its longi- 
tudinal compression energy. As a result, the phases of 
different segments are uncorrelated and the system again 
behaves as a collection of metallic rods. Polarizability of 
the rods generates a strongly anisotropic dielectric con- 
stant. Like for a single chain, in the ground state of the 
system, a finite bare density of states gs at zero energy 
originates from random background charges of the rods. 
Again this leads to the VRH at low temperatures. 

In the other case (Fig.^), the concentration of impuri- 
ties is small and chains are strongly interlocked. The elas- 
tic distortions are concentrated in small regions around 
individual impurities (see below). Away from impurities 
the crystal possess a good 3D order. The true long-range 
order is however absent because of the cumulative ef- 
fect of small elastic distortions in a large volume. The 
elastic displacement field u(r) of the electron crystal lat- 
tice away from impurities gradually varies in space. The 
length scale where its variation is of the order of a (vari- 
ation in cj) is of the order of unity) is referred to as Larkin 
length. The Larkin length is exponentially largefSi* 2 ^ ef- 
fectively infinite, because the Coulomb interactions make 
the crystal very rigid. This will be discussed in more de- 
tail in Sec. |V| 

The region near a typical impurity has the following 
structured On one side of the impurity the chain is com- 
pressed, which creates an excess negative charge; on the 
opposite side, it is stretched resulting in a positive charge 
of the same absolute value < e/2. The net charge of such 
a dipole is zero (together with the vacancy of charge — e 
and the acceptor of charge e). The characteristic length 
of the distorted region is of the same order as the length l s 



of the nonlinear topological excitation2a2Si2& of the pure 
system, the 27r-soliton. This is because the magnitude of 
local distortions in (f> are typically comparable in the two 
cases. The formal definition of l s is2£ 

l a = a±/y/a, (1) 

where a = Y±/Y x <C 1 is the dimensionless anisotropy 
parameter and Y x , Y± are the longitudinal and trans- 
verse elastic moduli of the electron crystal. The elastic- 
ity theory of the crystal will be discussed in more detail 
in Sec. Here we just mention that a varies from ma- 
terial to material, e.g., a ~ 10~ 4 in KCP and in CDW 
organics, a ~ 10~ 2 in blue bronzepi&ii 

At low temperatures all solitons are bound to acceptors 
with large binding energies comparable to the creation 
energy of a free soliton W ~ c 2 /kI s . In other words, 
there is a large energy gap for creating charge excitations. 
Nevertheless, as shown in Sec. E| finite gs at zero energy 
does exist in the case I ^> l s as well. It comes from 
rare clusters of several closely spaced impurities. Such 
clusters can be viewed as microscopic inclusions of the 
I -C l s phase (where gs is large). 

In all cases outlined above, gB is not yet the actual 
density of states of charge excitations. This is because 
the long-range Coulomb interaction of charges at distant 
sites is not included in the definition of gs- We denote 
the true density of states of charge excitations by g{e). 
Based on previous studies of other insulating systems, 
such as doped semiconductors^^ we expect that long- 
range interactions generate a Coulomb gap in g(s). This 
gap is soft, in the sense that 17(e) vanishes only at the 
Fermi level e = 0. Away from the Fermi level, g(e) in- 
creases in a power-law fashion until it saturates at the 
bare value at large enough e (cf. Ref. and Fig. 10.4 
therein) . Note that g is different from the thermodynam- 
ical density of states. The latter does not vanish despite 
Coulomb correlations, see Refs. 13 andE^. 

In macroscopically isotropic electron systems the func- 
tional form of the Coulomb gap depends on the number 
of dimensions. The density of states behaves as e 2 in 3D 
and as |e| in 2D. In all dimensions, however, this leads to 
the the Efros-Shklovskii (ES) law for the VRH conduc- 
tivity in isotropic doped semiconductors^* 2 ^ 

< 7 = ( T exp[-(T E s/T) 1 / 2 ], (2) 

where 00 is a prefactor, which has an algebraic In- 
dependence. Parameter Tes is given by 

T ES = Ce 2 M, (3) 

where £ is the (isotropic) decay length of localized elec- 
tron states, k is the (isotropic) dielectric constant of the 
semiconductor, and C is a numerical coefficient. (We 
measure temperature in energy units throughout this pa- 
per.) In lightly doped isotropic semiconductors n and £ 
are determined solely by material parameters (binding 
energy of impurity states, electron effective mass, band- 
structure, etc.). Therefore, Tes does not depend on the 
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impurity concentration and Eqs. J5J and Q are in this 
sense universal. 

In contrast, in this paper we show that in strongly 
anisotropic systems the Coulomb gap has, in general, a 
different functional form. Depending on I and other pa- 
rameters, it may be either universal or not [i.e., g(e) and 
Tes may contain factors related to impurity concentra- 
tion] . These results and their consequences for the VRH 
transport are presented in the next section. 



II. RESULTS 



of Ref. |3l| by two logarithmic factors. One of them orig- 
inates from logarithmic charging energy of metallic rod; 
the other — from the logarithm in the tunneling action 
s, see also Refs. ITil fl9l andl32l 

Once the logarithmic Coulomb gap is taken into ac- 
count, the T-dependence of the conductivity can still be 
written in the form of a ES law [Eq. (J5J] but Tes becomes 
a function of T, as follows: 
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We group our results according to the three cases (A, 
B, and C) outlined in the Introduction. 



A. Single chain 

In this case, studied in Sec. II I II the Coulomb inter- 
action is not screened. However, in ID the 1/x-decay 
of the Coulomb potential is on the borderline between 
the short and the long range interactions. Consequently, 
most physical quantities differ from their counterparts 
for the short-range (screened) interaction only by some 
logarithmic factors. For example, the density of states of 
charge excitations exhibits a logarithmic suppression, 30 
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ln(e 2 /^|e|) 



(4) 



In a first approximation, such a suppression can be disre- 
graded in the calculation of the VRH transport. Namely, 
one can assume that 17(e) = gg = const. In this approx- 
imation one arrives at the conventional Mott VRHf2*2& 
which in ID coincides with the ES law of Eq. J5J. Let 
us denote by T E ?^ the value of Tes that one obtains ne- 
glecting the Coulomb gap, then T^) ~ l/gs^x- Here 
£ x stands for the localization length that determines the 
asymptotic decay P oc exp(— 2x/^ x ) of the probability 
of tunneling of charge-e excitations over a large distance 
x. If the probability of tunneling between nearest rods 
is written in the form exp(— 2s), where s 3> 1, then 
tunelling paths with returns can be neglected. In this 
case £ x — l/s. Using the expression for s from Ref. IlH 
we obtain 
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In the last equation we absorbed numerical factors into 
the coefficient C\ ~ 1. 

A similar expression for T^ was obtained by Natter- 
mann et alm^ for the model of a disordered Luttinger 
liquid with short-range interactions and weak pinning. 
Our Eq. JSJ and © differ from the corresponding results 



Note that in the ID case, the standard derivation 3 , of 
the VRH law (J2J overlooks the role of very resistive hops 
in some exponentially rare places along the chain. A 
more careful approach shows 3 ^ 3 . that Eq. (J2J and its 
generalization through Eq. @ are valid only if the chain 
is sufficiently short. The quantitative criterion on the 
length of the chain can be obtained following the excellent 
discussion in Refs. |3j| and|32- This, however, goes beyond 
the scope of the present work. 



B. 3D systems with large impurity concentrations 

This case, formally defined by the inequality a± <C I <C 
l s is studied in Sec. IIVI It may be realized in strongly 
anisotropic CDW compounds such as KCP where the 
soliton length l s is large (10 2 a or so) and/or in sam- 
ples where a relatively high impurity concentration is 
created intentionally 34 35 so that I is small. Possible non- 
CDW realizations include arrays of relatively distant ID 

conductors, e.g., quantum wires, nanotubes, 4 — or poly- 

44 

mers^* 

As elaborated in Sec. ^ impurities divide the system 
into a collection of metallic rods. The finite 3D concen- 
tration of highly polarizable rods results in a large di- 
electric constant^ along the x-axis. The Coulomb inter- 
action is therefore strongly anisotropic but the Coulomb 
gap remains parabolic, g(e) oc e 2 , as in isotropic sys- 
tems. Tunnelling is anisotropic as well. The interchain 
tunnelling is accomplished by single-electron like excita- 
tions, which do not perturb charges on the intermediate 
chains along the tunneling path. In App.[B]we estimate 
the corresponding transverse localization length £j_ to be 
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ln(e 2 / natj 



(8) 



where t± is the interchain hopping matrix element in the 
tight-binding band-structure model. For the low-T VRH 
conductivity we again obtain the ES law with Tes given 
by 



Tes = C2 — 
kI 
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Here C2 is another numerical factor of the order of unity. 
We see that in both cases, A and B, the ES law looses its 
universality, because Tes depends on the impurity con- 
centration N through I = 1/Na\. For a single chain 
(case A) this dependence originates mainly from the de- 
pendence of the localization length £ x on I. In 3D (case 
B), the full effective dielectric constant and therefore, the 
density of states inside the Coulomb gap depend on N as 
well. In both cases, with decreasing N the temperature 
Tes decreases, which at a fixed temperature leads to an 
exponentially increasing conductivity. 

In doped semiconductors similar violations of the 
universality of Eq. @ are known to occur near the 
metal-insulator transition. In that case, however, Tes 
has an opposite dependence on N. In particular, 
Tes vanishes when N grows and reaches the criti- 
cal concentration^ Similarly, previous theories of the 
VRH transport in strongly anisotropic systems dealt 
with gapped, semiconductor-like materials (commensu- 
rate CDW) where impurities provided carriers f^LS so 
that the conductivity was found to grow with the impu- 
rity concentration. In contrast, our work is devoted to 
systems, which are metallic (sliding) in the absence of 
impurities. Therefore, decrease of Tes with decreasing 
TV seems naturali^S 

At high temperatures, the conductivity is due to the 
nearest-neighbor hopping. Its T-dependence is of acti- 
vated type, 
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a = g a cxp(-T j4 /T), I < l s , 
with the activation energy 



A kI 



(10) 



(11) 



and the prefactor a a proportional to the probability of 
tunneling between adjacent rods. 



C. 3D systems with small impurity concentration 

As impurity concentration decreases and I becomes 
larger than l s , a number of dramatic changes appear in 
all key quantities, such as the density of states, the lo- 
calization length, and the effective dielectric constant. 
For example, as we discuss in Sec. [V] the dielectric con- 
stant starts to increase exponentially with I because the 
polarizability of the crystal with / > l s becomes lim- 
ited not by the length I of individual chain segments 
but by the exponentially large length of Larkin domains. 
The soaring dielectric constant causes a rapid drop of 
the ES parameter Tes- In turn, this causes a collapse 
of the low-temperature resistivity in a narrow interval 
Is i$ I l a \x\.(yV/T) (see the descending branch of the 
curve in Fig. 0). Until this point, the notion that our 
system is opposite to the conventional semiconductors, 
so that purer samples have higher conductivities, seem 
to be working. 
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FIG. 2: Logarithm of the resistivity as a function of the aver- 
age inter-impurity distance I = 1 /Na± at a fixed temperature 
T -C W. The ES, 2/5, Mott, and activation laws succeed each 
other with growing I. 



Once I exceeds l s ln(W/T), the Larkin length L x can 
be treated as effectively infinite. The VRH now involves 
hops between low-energy states separated by distances 
shorter than L x . On such distances, the dispersion of 
the dielectric function becomes important. Each pair of 
low-energy charge excitations localized on their respec- 
tive impurity clusters interacts via a strongly anisotropic 
electrostatic potential, which is not exponentially small 
only if the vector that connects the two charges is nearly 
parallel to the chain direction. Such an unusual inter- 
action leads to a Coulomb gap that is linear in energy 
and independent of N, unlike the previous case (case B, 
I <C l s ), where the Coulomb gap is quadratic and N- 
dependent. Another difference from the case B is that 
the localization length £ x for the tunnelling in the chain 
direction is also independent of N, 



Is 



(12) 



see Sec.lVland App.[Bj This leads to a novel 2/5-law for 
the VRH conductivity, 



< 7 = ( roexp[-(T 1 /T) 2 / 5 ], 



where parameter Ti, given by 



T 1 = C l 



2 1/4 
e TV aj_ 

kI s £_l 



(13) 



(14) 



does not depend on I and is, in this sense, universal. 
Here C3 is yet another numerical coefficient of the order 
of unity. The 2/5-law shows up as an intermediate resis- 
tivity plateau in Fig. This universal law for quasi- ID 
systems with / 3> Z s is an analog of the universal ES law 
in isotropic systems. 

We show in Sec.0that the Coulomb gap affects mainly 
a finite energy interval |e| < A, where A cx gs can be 
called the Coulomb gap width. At larger energies, the 
density of states of charge excitations coincides with the 
bare one, g{e) ~ gs- Since gg is generated by impu- 
rity clusters whose concentration diminishes with grow- 
ing I, both gB and A decrease with /. Eventually, the 
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Mott 



FIG. 3: Summary diagram for the transport regimes in a 3D 
system. Domains of validity of ES [Eq. ©], Mott [Eq. ITSt l. 
activated [Eqs. I|1O0 and l|16|l ] and 2/5 [Eq. I|13|l ] laws are 
shown. 



Coulomb gap becomes more narrow than the range of 
energies around the Fermi level responsible for hopping 
at given T. At this point the Coulomb gap can be ne- 
glected and the 2/5-law is replaced by the conventional 
Mott law, 



<t = <t exp[-(T M /T) 



l/4l 



(15) 



where Tm = Ca/ gB£,x£,±, Ca ~ 1. As I increases further, 
T being fixed, a decreases because of diminishing gB- 
This gives rise to the ascending branch in Fig. El At such 
I, the electron crystal behaves similar to a gapped insula- 
tor where a lower impurity concentration corresponds to 
a lower carrier density, and thus, to a higher resistivity. 

As I continues to grow, at some point the Mott VRH 
crosses over to the nearest-neigbor hopping and shortly 
after it becomes smaller than the conductivity due to 
thermally activated free solitons, 



a = a A csq>(-W/T). 



(16) 



At even larger I, a ceases to depend on I, and so the 
impurity concentration, see Fig. [21 Note that the acti- 



vation energies W ~ e 2 / ' nl s [Eq. (TSJ] and Ea 



1 /kI 



[Eq. I|lll) ] in the cases A and B, respectively, smoothly 
match at I *~ L. 



D. Summary of the regimes 

The rich behavior of the conductivity as a function of 
I and T is summarized in the form of a regime diagram 
in Fig. |31 The novel 2/5-law applies in a broad range of 
I and T between the ES and the Mott laws. 

A convenient way to keep track of all the VRH expo- 
nents derived in this paper is provided by Eq. (|17f) below. 
We would like to present it in a somewhat more general 
form, motivated by the following physical reasoning. 

The diagram of Fig. |3] is obtained under the assump- 
tion that the tunneling in the transverse direction is not 
negligible is not too small), so that the VRH has 
a 3D character. However, if the conducting chains are 



TABLE I: The exponents A of VRH conductivity [Eq. 1X71 1 
in the cases of 3D, 2D and ID tunneling and a power-law de- 
pendent density of states g(e) that arises due to 3D Coulomb 
interactions. ID tunneling corresponds to £x — » 0. 
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relatively distant from each other either along g or in- 
direction or both, this condition may be violated. Ex- 
amples of such systems are artificial arrays of quantum 
wiresi and carbon nanotubes^ In those systems, the 3D 
hopping is pushed to very lower temperatures, while at 
intermediate T the hopping can be either one or two- 
dimensional. Generalizing the standard derivation of the 
VRH transport^* to the <i-dimensional hopping and a 
power-law density of states g(e) oc e M , one obtains the 
conductivity in the form 



cr = Co ex P 



- (TvRh/IT 



A = 



M+ 1 



1 



(17) 



For d = 3 one recovers all the regimes discussed prior in 
this Section (Mott, 2/5, and ES laws) by setting /j, succes- 
sively to 0, 1, and 2, according to the physical situation. 
For the sake of completeness, the exponents for other d's 
in the same situations are summarized in Table [TJ In- 
clusion of all such regimes would transform Fig. into a 
more complicated diagram, but would not change its gen- 
eral structure, so it will not be shown here or discussed 
further below. 



III. ID SYSTEM 

In the case of a single chain of electrons on a uni- 
form positive background (Fig. [TJi) impurities divide the 
ID electron crystal in separate pieces, which behave as 
metallic rods. The rod lengths x are distributed ran- 
domly around the average value I. Therefore, the back- 
ground charge of a given rod, Q — —ex/a, is a random 
number. It can be written as Q = —e(n + u), where n 
is an integer and v is a number uniformly distributed in 
the interval -1/2 < v < 1/2. 

In the ground state of the system each rod contains 
an integer number n r of electrons, so that the rod has 
the total charge of q = e(n r — n — v). To find n we 
use the fact that the Coulomb self-energy of the rod is 
equal to q 2 /2C r where C r = Kx/[2ln(x/a)] is the capac- 
itance of this rod. On the other hand, the interaction 
of different rods does not contain the large logarithm 
ln(x/a), and can be neglected in the first approximation. 
Thus, the minimization of the total energy of the sys- 
tem amounts to minimizing the self-energies of the rods. 
One can show then that, in the ground state, n r = n, 
so that the charges of the rods are uniformly distributed 
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in the interval — e/2 < q < e/2. Indeed, if this is not 
true and the charge of the given rod is q > e/2, then, by 
charge conjugation symmetry, there should exist another 
rod with the same length and the opposite charge —q. 
Transferring an electron from the first rod to the sec- 
ond one we lower the total energy and make the absolute 
values of both charges smaller than e/2. 

Below we use q~—ev and call rods with < v < 1/2 
empty and with — 1/2 < v < occupied. We define the 
energy of an empty state, e(x, v), as the minimum work 
necessary to bring to it an electron from a distant pure 
ID electron crystal with the same average linear density 
of electrons 



e(x, v) 



In 



KX 
„2 



(I) 



-ln(-) (1-2*), 
kx \a/ 



(18) 



This energy is positive and vanishes only at v = 1/2. On 
the other hand, the energy of an occupied state is defined 
as minus the maximum work necessary to extract elec- 
tron from this rod to the same distant pure ID electron 
crystal. In this case the final result is identical to Eq. I|18|) 
with v — > — v. Apparently, states with \v\ = 1/2 are ex- 
actly at the Fermi level, which we take as the energy 
reference point (it coincides with the electron chemical 
potential of a pure ID crystal). The low-energy states 
relevant to VRH transport have \v\ — 1/2 1. 

Now we can calculate the density of such states. Tak- 
ing into account the fact that this rod length x is dis- 
tributed according to Poisson statistics, the disorder- 
averaged bare density of states can be written as 



9b(E) 



1/2 oo 

f dv I ^ exp (-^js(E-e(x,v)). (19) 



With a logarithmic accuracy, we can replace hi(x/a) by 
\n{l/a). Then we use Eqs. (JTHJl and (JTHJ) to find 



9b(s) 



1 

2eal 



1 



1 



so 



exp 



(-?) 



(20) 



where e — e 2 \ti(1 / a) / nl . We warn the reader that one 
should not attribute much significance to the predictions 
of the above formula in the region of high energies, e > Eo, 
where excitations with charges larger than e will also 
contribute to various physical processes. On the other 
hand, close to the Fermi level, for e <C £q, only charge-e 
excitations are important, in which case Eq. I|20|) is fully 
adequate, while gs is nearly constant, 
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2e 2 \n(l/a)' 



(21) 



For the calculation of the low-T transport, only this con- 
stant value is needed. 

As mentioned in Sec.^ in ID the l/x-Coulomb inter- 
action creates only marginal effects on the conductivity. 



If we neglect them, in the first approximation, then stan- 
dard Mott's argument*? leads to the VRH that obeys the 
T 1/2 -law, Eq. ©, with T ES = T M = C x /g B i x , where 
£ K is localization length for tunnelling between distant 
rods. The value of £ x is obtained from the following con- 
siderations. Tunneling through an impurity that sepa- 
rates two adjacent rods can be viewed as a process in 
imaginary time that consists of the following sequence of 
events^ A unit charge assembles into a compact soliton 
just before the impurity in one rod, tunnels through the 
impurity, and finally spreads uniformly over the other 
rod. We assume that the chain is not screened by exter- 
nal metallic gates. Then the tunneling probability can 
be written in the form exp(— 2s), where the s is the di- 
mensionless action^ s ~ r X J 2 ln 3 / 2 (Z/a). Since s 1, 
tunneling paths with returns can be neglected. There- 
fore, for electron tunnelling over distances x 3> I the 
action s should be multiplied by the number of impuri- 
ties passed. The average number of such impurities is 
equal to x/l, which yields the total tunnelling probabil- 
ity P cx exp(— 2x/£ x ), with £ x given by Eq. (JSJ. With 
a logarithmic accuracy, the effect of the Coulomb gap 
is to replace qb by g[e) evaluated at the characteristic 

hopping energy e = \J T^T, where is defined by 
Eq. ©. The latter result follows from Eq. (J2J and JSJl. 
The final expression for the parameter Tgg is given by 
Eq. J2J|. We briefly note that in a long enough chain 
clusters of atypically densely spaced impurities may ex- 
ist. Tunneling through such segments would cost a higher 
tunneling action and therefore the overall conductivity 
would be suppressed. In this paper we assume that the 
chain is sufficiently short so that these rare clusters can 
be neglected, see a comment after Eq. {7J). 

Formulas © and JJJ indicate that Tes goes down as 
impurity concentration N = 1/1 decreases. This provides 
a gradual crossover to the metallic behavior in a pure ID 
system. 



IV. 3D SYSTEM WITH A LARGE IMPURITY 
CONCENTRATION 



In this section we consider a quasi-lD system made 
of parallel chains which form a periodic array in the 
transverse directions. The chains are pinned by impurity 
centers, which divide them into metallic rods of average 
length / = 1/Naj_, where a\ is the cross-sectional area 
per chain. The finite-size electron crystals in each rod 
are either compressed or stretched to accommodate an 
integral number of electrons, as in the case of a single 
chain. In this section we assume that for a typical rod, 
the energy of its longitudinal deformation is smaller than 
the energy of its transverse coupling to rods on the neigh- 
boring chains. This can be the case when the periodic 
potential created by the neighboring chains is diminished 
because a± is larger than a, and/or when the impurity 
concentration is large enough. Formally, the inequality 
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aj_ <C I <C l s needs to be satisfied, where l s is the soli- 
ton length [Eq. 0]. Since metallic rods now completely 
fill the 3D space (see Fig. EJ)), they modify the dielec- 
tric constant of the system. As in the interrupted-strand 
models the dielectric constant is anisotropic. Along the 
chain direction it has the value of 



k x = K[l + C 5 (l/r D )' 



(22) 



where rjj ~ a± is the screening length, see Apps.IaI and 
Cs ~ 1 is a numerical constant. 39 Transverse dielectric 



constants are unaffected, K y 



k. At large I, k x 



greatly exceeds K, which leads to an anisotropic Coulomb 
interaction in the form 4 ^ 



U(r) 



Ky/x 2 + (k x /k)t 2 ± ' 



(23) 



where r\ 



y 2 + z 2 . In spite of the large dielectric con- 
stant, the Coulomb interaction is long-range and thus 
creates a Coulomb gap. Using the standard ES argu- 
ments^ the following density of states of charge excita- 
tions is obtained 

gi£ )=*^, (24) 

It differs from the conventional formula for an isotropic 
medium only by the presence of k x instead of k. In or- 
der to calculate the VRH conductivity we still have to 
discuss the tunnelling probability. It is important that 
Eq. (|23|) holds only at x 3> l- The interaction between 
charge fluctuations on the same rod is short-range due to 
screening by neighboring chains. 45 Tunnelling along the 
x-axis takes place similarly to the case of a single chain, 
but screening of the Coulomb interaction leads to smaller 
action of the order of s ~ y/rJlnQ/a). Therefore, for the 
localization length in the x-direction we get 



I 



\n{l/a)' 



(25) 



The tunneling in the y- and z-directions is accomplished 
by single-electron like excitations. The probability of 
tunneling decays exponentially, P cx exp(— 2rj_/£x) ; a ^ 
large transverse distances r±. Here £,± is the transverse 
localization length given by Eq. JSJ and discussed in more 
detail in App. [B] We assume that £j_ is not vanishingly 
small compared to £ x , in which case the VRH has a 3D 
character and can be calculated by the percolation ap- 
proach (see Ref. |jand references therein) 4i This calcu- 
lation differs from the isotropic case by the replacement 
of the isotropic dielectric constant n and the isotropic lo- 
calization length £ by their geometric averages over the 
three spatial directions: 



Tes — C*3 



f 



(26) 



With the help of Eqs. l|2*2*)l and l|2*5|l the expression for 
Tes reduces to Eq. where the screening length of the 
electron crystal has been taken as rn ~ a± and numerical 
coefficients absorbed in C3. 



V. 3D SYSTEM WITH A SMALL IMPURITY 
CONCENTRATION 

In this section we study the crystal pinned by impu- 
rities with a low concentration iV so that the condition 
I ^> l s is satisfied. As discussed in Sec. [Hall key quantities 
— density of states, localization length, the screening of 
Coulomb interactions — undergo dramatic changes com- 
pared to the case of high impurity concentration. Wc 
address such changes in the three separate subsections 
below. 



A. Pinning of a quasi-lD crystal by strong dilute 
impurities 

In this subsection we study the ground state structure 
and screening properties of the crystal with low impurity 
concentration. 

We start by reviewing the physical meaning of l s given 
by Eq. 0J, in which a = Y±/Y x is the anisotropy pa- 
rameter, and Y x , Y± parametrize the energy of an elastic 
distortion of the crystal, 

E ° l = ^f d 3 r[Y x (d x u) 2 + Y ± (X7 ± u) 2 }. (27) 

As shown in App. at large r s the longitudinal elastic 
modulus Y x is dominated by Coulomb effects, 



v 2/22 

Y x ~ e / Ka a_L- 



(28) 



The transverse modulus Y± can be substantially smaller 
than Y x even when the ratio a±/a is only modestly large. 
For Coulomb interaction, lj_ cx exp(— 2na±/a). 

The energy E el in Eq. (|27|l is essentially the short- 
range part of the Coulomb energy. The total energy also 
includes the long-range Coulomb part (see App. 01 and 
the pinning part (see App.lcland below). Strictly speak- 
ing, Eq. (|27|l is valid only for small gradients of the elas- 
tic displacement field u(r); however, it can be used for 
order-of-magnitude estimates down to microscopic scales 
r±_ ~ a±_ and x ~ l s . In this manner, one can derive 
formula by minimizing E cl under the condition that 
u changes from to a over a segment of leng th l s on 
a single chain. For more details, see Refs. |23, 123, and 
App.|H 

The inequality I 3> l s imposes the upper limit on the 
impurity density: 



N < y/a/c 



(29) 



Below we show that at such N the ground state of the 
electron crystal is determined by an interplay of individ- 
ual and collective pinning4& 

Without impurities the crystal would have a perfect 
periodicity and long-range 3D order. Impurities cause 
elastic distortions of the lattice. The strongest distor- 
tions, of dipolar type, are localized in the vicinity of 
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impurities, see Sec. H] and Fig. QJ;. Such dipoles have a 
characteristic size l s (same as free solitons) , are well sep- 
arated from each other, and occupy only a small fraction 
of the space. Their creation is advantageous because the 
associated energy cost (elastic plus Coulomb) is of the 
order of W ~ e 2 / nl s per impurity, whereas the energy 
gain is a much larger electron-impurity interaction en- 
ergy — e 2 / 1 na. This is the essence of individual (strong) 
pinning phenomenon, which provides the dominant part 
of the pinning energy density. The collective (weak) pin- 
ning results from interaction between the dipoles. Let us 
demonstrate that such interaction cannot be neglected at 
sufficiently large length scales. By solving the elasticity 
theory equations (generalized to include the Coulomb in- 
teractions), it can be shown^l that a dipolar distortion 
centered at a point [Xi, rx») has long-range tails that de- 
cay rather slowly with distance, u ~ Aiaro/\/a\x — Xi\. 
This displacement is confined mainly within a paraboloid 
\r± — r±i\ 2 < y/ar£>\x — Xi\. Note that the segment 
< x < x min = ax (J/ ' \farr,) 1 !' 1 of the paraboloid 



'arnm 



(30) 



contains on average one impurity. Parameter ~ a±, 
which we already encountered in Sec. llVI has the meaning 
of the screening length. It is related to the longitudinal 
elastic modulus Y T as follows: 



' D 



Aire- 



: a 2 a 4 L Y x , 



(31) 



see App. [XJ 

If, in the first approximation, we choose to neglect the 
interaction among the dipoles, then we should simply 
add their far elastic fields treating the amplitudes — 1 < 
Ai < 1 as random variables. We immediately discover 
the logarithmic growth of u with distance, 



x *J~a.r D x 



([u(x,0) - u(0,0)} 2 ) 



N / dx' 



dr 2 ' ^ 



where r"" n = (v^rD^min) 1 / 2 ■ 

Alternative derivation of Eqs. Q32J1 and (|33|l based on 
energy estimates is given in App.|0 It elucidates that the 
slow logarithmic growth of (Am 2 ) is rooted in the impor- 
tant role of long-range Coulomb interaction in the elastic 
response of a quasi- ID crystal. An isotropic electronic 
crystal adjusts to pinning centers primarily by means of 
shear deformationaSi that do not cost much Coulomb en- 
ergy. As a result, in the isotropic crystal Am grows alge- 
braically with distance. In contrast, in quasi-lD crystals 
and CDW, where elastic dispacement is a scalar (elec- 
trons move only along the chains), no separate shear 
deformations exist. The build-up of the Coulomb en- 
ergy that accompanies longitudinal compressions trans- 
lates into an exceptionally large rigidity of the electron 
lattice and exponentially large L x and L± . 

At distances exceeding the Larkin lengths the dipoles 
can no longer be treated as independent. Indeed, the en- 
ergy cost E s of a given dipole is determined by the mini- 
mal distance by which the crystal has to distort to align 
an electron with the impurity position. Therefore, just 
like the energy of a rod in the previous sections, E s has 
a periodic dependence on v — {{xi — u)/a}, where {. . .} 
denotes the fractional part. E s vanishes at v — and 
reaches a maximum value of~Watz/ = ±l/2. There- 
fore, as soon as the cumulative effect of other dipoles 
attempts to elevate \v\ above 1/2, a 2-7r-phase slip should 
occur to adjust E s (u) to a lower value. The overall ef- 
fect of such adjustments is to enhance the pinning en- 
ergy. This additional energy gain can be viewed as the 
collective ■pinning effect. Using standard arguments (see 
App. [UJ, we relate the extra pinning energy density to 
the Larkin length 



£ 



pin 



-W/a\L x ~ -Y^a/L^f 



(34) 



We will now use this result to estimate the asymptotic 
value 



efe -> 0, q± = 0) 



(35) 



V\ ln ( x 

C*6 J \ l J \X m in 



(32) 



where the bar over u indicates that we refer to the value 
of u away from the immediate vicinity of a dipole and 
(...) stands for disorder averaging. 

The logarithmic growth of the elastic displacements 
was previously derived for the model of weak pinning 
centers in early worksSi^ on the subject. In those cal- 
culations the numerical coefficient Cq is large and is in- 
versely proportional to the impurity strength. In our case 
Cq ~ 1. Apart from that, Eq. 1)32[) demonstrates that 
the case of strong pinning centers is essentially similar. 
Therefore, as customary for weak pinning models we de- 
fine the longitudinal L x and transverse L± Larkin lengths 
as the lengthscales where (Am 2 ) = ([it(r) — u(0)] 2 ) ~ a 2 . 
From Eqs. 1(50)1 and we obtain 



of the longitudinal component of the dielectric function. 

Without impurities, the dielectric function has the 
form 



*(q) 



• D 



q 2 q 2 + aq^ 



(36) 



see App. El an d, e.g., Ref. 0- Based on previous 
wor j i io J 22 J 48 we assume that a reasonable description of 
dielectric screening in a system with impurities is ob- 
tained if we replace the random distrubution of pinning 
centers by a commensurate pinning with the same £ p j n . 
In this case we can use the concept of the dielectric func- 
tion even for the disordered system. To derive the mod- 
ified expression for the dielectric function, one can add 



the term £ p i n (u/a) to the right-hand side of Eq. (|A2I) in 
App. El an d repeat the steps outlined thence. It is easy 
to see that the net effect of (commensurate) pinning is 
' enp(Cel/2l s ), (33) to augment the combination q 2 + aq\ (proportional to 
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the elastic resoring force) by the term — £ p i n /Y±a 2 , which 
comes from the additional restoring force due to impuri- 
ties. In this manner we obtain 



e(q) 



q 2 q 2 x 



td 



LI 2 ) 



= —;= exp 6 6 - 7 



(37) 
(38) 



However, for our purposes a cruder approximation will 
be sufficient. Namely, we can assume that at distances 
shorter than the Larkin length, the system screens as 
though it is free of impurities, Eq. I|36(l : at distances 
larger than the Larkin length, the dielectric function is 
replaced by a constant, Eq. QH8JI. In the following subsec- 
tion we will use Eqs. I|3t)|) and l|38|l to derive the functional 
form of the Coulomb gap in the regime I l s . 



Bare density of states and the Coulomb gap at 
low impurity concentration 



b) 



o) 



2 M+l M+2 



(i + f)e 



(| + f)e 



(i + |)e 



2 M+l M+2 



a-S)e 



2 M+l M+2 



a-S)e 



2M+2 



(i + f)e 



2M+2 



H + f)e 
2M+2 



(| + 



FIG. 4: An example of an impurity cluster that provides 
low-energy charge excitations at I 2> l s - (a) Distribution of 
charges in the ground state. For simplicity, we chose 5l = 
8/2. The cluster is electrically neutral, (b) Competing neutral 
state, (c) Excited state obtained by increasing the charge of 
the central rod by one unit. 



In order to describe the low-T transport at low impu- 
rity concentration, I 3> l s , we need to determine the ori- 
gin of low-energy charge excitations in this regime. This 
poses a conceptual problem. Indeed, such excitations do 
not exist in the bulk (away from impurities) where the 
creation energy of charge-e excitations is bounded from 
below by the energy W of a 27r-soliton. At first glance, 
the impurities do not help either. As mentioned in SecQJ 
near isolated impurities there is an energy gap for charge 
excitations, which is not much smaller than W. This is 
because a single impurity appreciably disturbs the crys- 
tal only within the region of length l s . The disturbance 
is electrically neutral (dipolar) in the ground state. Cre- 
ation of a charge-e excitation near such an impurity re- 
quires an energy of the order of e 2 / nl s ~ W. Let us now 
show that charge excitations of arbitrary low energies 
nevertheless exist. They come from impurity clusters. 
Each cluster is a group of a few impurities spaced by dis- 
tances of the order of l s or smaller. (It can be viewed as 
a microscopic inclusion of the I < l s phase.) Below we 
demonstrate that the clusters provide the bare density 
of states at zero energy, which decreases with I no faster 
than a power-law, 

9b= ^[j) ' (39) 

where exponent (3 is of the order of unity and is indepen- 
dent of The calculation of (3 has to be done numerically, 
which we leave for future work. 

To prove that Eq. gives the lower bound on gg, 
consider the configuration of 2M + 2 impurities shown in 
Fig. 0] The impurities define a cluster of 2M + 1 short 
rods, each of approximately the same length c <ti l s - The 
cluster is flanked by two semi-infinite segments at the 
ends. We assume that M is sufficiently large so that 



L = Mc is much greater than l s . Suppose that the length 
of the central rod in units of a is close to a half-integer, so 
that the charge of this rod is restricted to the set of values 
q = (— 1/2 — £+n)e, where n is an integer and < 5 <C 1. 
For the low energy states we only need to consider two 
possibilities, q — {— 1/2 — 5)e and (1/2 — 5)e. The lengths 
of the other short rods in our construction are chosen 
to be close to integer multiples of a. Then those rods 
can be considered charge neutral. Finally, we assume 
that the position of the leftmost impurity restricts the 
charge of the left semi- infinite cluster to 1/4 + Sl + til, 
where \8l\ <C 1 and nj, is another integer. Under these 
assumptions, the charge of the right semi- infinite segment 
is fixed to the set 1/4 + 5 — 5l + tir. 

One possible candidate for the ground state is the 
charge configuration shown in Fig. It is overall neu- 
tral, has the charge of the central rod equal to q = 
(—1/2 — S)e, and the total energy 

Ea —ni\r2) + 2rAv + — 2o — + 2Wl/4 

- -|£^(V2 + *) 2 + 2IFV 4 - (40) 

Here we made use of the conditions L > l„ i5 C 1. 
W1/4 denotes the self-energies of the semi- infinite seg- 
ments (which contain "1/4-solitons" ), and C r ~ kc is 
the self-capacitance of the central rod. If the configura- 
tion shown in Fig. is indeed a ground state, it has an 
unusual feature that the charge of the central rod exceeds 
1/2 by the absolute value, contrary to the arguments in 
Sec. IIIII To establish the conditions for the cluster to 
have such a ground state, we need to compare E a with 
energies of other possible states. Of those, some are neu- 
tral and some are charged. One competing neutral state 
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is shown in Fig. ^J). It has energy 

Eb = ~^h + ^c r il/2 ~ 6)2 + Wl/4 + W3/i ' (41) 

where W3/4 is the energy of the "3/4-soliton" on the right 
of the cluster in Fig.^J). For the energy difference we have 



Eh — E„ 



1 e z 



Wo 



C r - 8 Lk '^/4 "^1/4, 

so that the state shown in Fig. wins if 



S< 



(W 3/4 - W 1/4 ) 



8 Lk 



(42) 



(43) 



Since W 3 / 4 — Wi 



1/4 



W ~ e 2 /nl s and L » Z s , the right- 
hand side of the inequality l|4^|) is positive, and so 5 > 
that satisfy such an inequality do exist. The only other 
vaible competing neutral state is similar to that shown 
in Fig. except the '3/4-soliton" is formed on the left 
semi- infinite segment. The energy of that state is also 
Eh, and so it does not lead to any further restrictions on 
S. 

Now let us examine the charged states. There is only 
one viable competitor, with q — (1/2 — 5)e, as shown in 
Fig. If the energy difference e = E c — E a is positive, 
then Fig. 0k represents the true ground state and e gives 
the creation energy of the charge-e excitation. If e is neg- 
ative, then the ground state is as shown in Fig. Qt, while 
the charge and the energy of the lowest energy charge 
excitation are equal to — e and — g, respectively. An ele- 
mentary calculation yields 



a 



■5 



e 

2kL' 



(44) 



from which we conclude that it is possible to obtain arbi- 
trary small e of both signs by tuning S sufficiently close 
to C r /2LK r~j c/L, without violating the inequality (|43|l . 
This proves that clusters are able to produce a finite 
<7b(0). In the ground state some of these clusters are 
neutral ( "empty" ) and some are charged ( "occupied" ) . 

Now let us try to estimate gs due to clusters. In the 
above argument we required strong inequalities c« I s 
and L 3> l s to prove the existence of a nonzero gs with 
mathematical rigor. Physically, it seems reasonable that 
such inequalities can be softened to c < l s and L > l s , in 
which case 6 < 1 and M ~ 1, i.e., only a few (of the order 
of unity) impurities are needed. It is also clear that the 
rods to the left and to the right of the central one need 
not be exactly neutral for the argument to go through. 
Then the probability of forming the desired cluster is 
comparable to the probability of finding 2M + 2 (a few) 
impurities on the same chain with nearest-neighbor sepa- 
ration less than l s . Assuming that impurity positions are 
totally random and independent, we obtain the estimate 
of Eq. Qg). 

Now let us calculate the form of the Coulomb gap in 
the actual density of states 17(e). At exponentially small 



energies, e <C e* = e 2 /nL x , the Coulomb gap is deter- 
mined by the interactions on distances exceeding the size 
of the Larkin domain. At such distances the interaction 
has the form (|23|) , leading to the usual parabolic Coulomb 
gap given by Eqs. lj2*4l and (J3SJ. We put these equations 
side by side below for the ease of reading: 



9(e) = ~ 



3 K K x 2 

7r e u 

a 2 a 2 1 



kL x y/ar D 



exp ( C 6 — 



(45) 
(46) 
(47) 



To ascertain the region of validity of Eq. I|45[l one needs 
to make sure that g(s) does not exceed the bare density 
of states gs- It is easy to see that this is the case here. At 
the largest energy e ~ [Eq. (|4^1) ]. g(e) is exponentially 
small, g oc exp(— Cel/l g ), and so it is indeed much smaller 
than gs ~ (l s /l)^ +2 - In this sense gs is large enough to 
ensure the validity of the parabolic law over the full range 
of e specified in Eq. i|43|l . 

At e larger than e», the Coulomb gap is governed 
by interactions within the volume of a Larkin domain 
and the dispersion of the dielectric function e(q) be- 
comes important. The interaction potential is defined 
by U(q) = Aire 2 /e(q)q 2 and Eq. (|3TJl in the g-space. In 
real space, it is given bySl 



U(v) 



2k\x 



exp 



A\fa ro\x\ 



si < ri < min{ v / a|a;|, L±}, 



(48) 

< x < L x . (49) 



The potential U(r) is appreciable only within the 
paraboloid defined by Eq. (j20J), where it behaves as 
U(r) ~ e 2 /2n\x\. A more precise statement is that the 
surface of a constant U is a paraboloid-like region 



4v / aru|a;| In 



2Uk\x\ 



(50) 



For U ^> e*, this surface belongs to the domain l|49|) 
where Eq. holds. 

To calculate the functional form of the Coulomb gap we 
use the self-consistent mean-field approximation due to 
Efros, 42 according to which g is the solution the integral 
equation 



g{s) = 9b exp 



w 



~\ I <l'<li.:')Vi: :') 



(51) 



where V(U) is the volume enclosed by the constant-C/ 
surface (|5Up. 



V(U) 



4 n 2 U 2 ' 



(52) 
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FIG. 5: The density of states of charge excitations in a 3D 
system with I 3> h- The parabolic Coulomb gap at low ener- 
gies is succeeded by the linear Coulomb gap, then by the bare 
density of states gs created by impurity clusters. 



Equation i|51|) follows from the requirement that the 
ground state must be stable againts a transfer of a 
unit charge from an occupied state with energy — e' to 
an empty state with energy e. Such a stability cri- 
terion*^ can be expressed by means of the inequality 
e + e' — U(r) > 0, where r is the vector that connects the 
two sites. Thus, the integral on the right-hand side of 
Eq. (|51|l counts the pairs of states that would violate the 
stability criterion if positioned randomly. The Coulomb 
gap can be viewed as the reduced statistical weight of 
stable ground-state charge distributions with respect to 
that of the totally uncorrelated ones. 

The solution of Eq. I|51[) has the asymptotic form 



9(e) 



8 



TT^/a roe 4 



r* < \e\ < A. 



(53) 



We see that the unusual interaction potential of Eq. I)48|) 
leads to a nonstandard Coulomb gap, which is linear in 
a 3D-system. The weaker (linear instead of the standard 
quadratic) suppression of g is due to a metallic screening 
of the Coulomb potential in the major fraction of solid 
angle. The energy scale A in Eq. I|53[l is defined by the 
relation g(A) ~ gB- Using Eq. I|39|) . we can estimate A 
as follows: 



-e gB 



W 



0+2 



(54) 



At e 3> A, the solution of Eq. I|51|> approaches the bare 
density of states, g ~ gs, and so A has the meaning 
of the Coulomb gap width. On the lower-energy side, 
at e ~ £*, the linear Coulomb dependence of Eq. 
smoothly matches the quadratic dependence of Eq. (|45|l . 
All such dependencies are summarized in Fig. [5] 



VRH transport in a system with small impurity 
concentration 



The VRH transport at Z 3> Z s involves quantum tunnel- 
ing of charge excitations between rare impurity clusters. 
In App.[B]we advance arguments that the charge tunnels 



along the chains in the form of 27r-solitons and derive the 
following estimate for the corresponding longitudinal lo- 
calization length: 



< .A. 

Sa; , 



Is i / Is 

-In - 

/ \ a 



I > L 



(55) 



At I ^ l s this formula goes over to Eq. I|12fl . while at 
I ~ l s it smoothly matches with Eq. I|25|) . As for the 
interchain tunneling, it is still accomplished by single- 
electron like excitations and the corresponding localiza- 
tion length £j_ is still given by Eq. (jHJ), see Add. IbI We 
will assume that £j_/a_i_ is not vanishingly small, in which 
case the VRH has the 3D character^ 

The density of states and localization lengths is all the 
information we need to calculate the VRH conductivity. 
A more complicated dependence of the density of states 
g(e) on energy in the case at hand, I > l s , brings about 
a larger variety of possible transport regimes. 

At the lowest temperatures we still have the ES law 
with parameter Tes given by Eq. \2ty . Due to the 
exponential growth of the longitudinal dielectric con- 
stant k x with I [Eq. I|38|) ]. Tes decreases exponentially, 
Tes oc exp(— CqI/31 s ). This dependence entails a precipi- 
tous drop of — In a as a function of I, as soon as I exceeds 
l s , see Fig. [5] Such an enhancement of the conductiv- 
ity is due to progressively more efficient screening of the 
long-range Coulomb interactions (steeply increasing k x ), 
which enhances the density of states inside the Coulomb 
gap, see Eq. tlSl . 

The ES law holds until the range of e's that con- 
tribute to the VRH transport^ |e| < (TT ES ) 1/2 , fits in- 
side the parabolic part of the Coulomb gap, |e| < e*. For 
a fixed T < W this gives the condition I < l s ln(W/T). 

At larger I, the unusual linear Coulomb gap l(53*|) leads 
to the novel 2/5-law for the VRH transport [Eq. lfT3))]. 
which we reproduce below for convenience: 



a = ( roexp[-(T 1 /T) 2 / 5 ]. 



(56) 



As emphasized in Sec. \n\ parameter Ti is impurity- 
independent and is, in this sense, universal, see Eq. i|14[l . 
This behavior leads to the intermediate plateau at the 
graph in Fig. [2] 

The range of energies that contributes to the VRH 
in the 2/5-law regime is given by |e| < (Ti/T) 2 / 5 T. 
At I ~ l s (W/T)i, where 7 = 3/[5(/3 + 2)] < 0.3, this 
range becomes broader than the Coulomb gap width A 
[Eq. I|54|) ]. At such and larger I, the Coulomb gap can be 
neglected, and the VRH begins to follow the usual Mott 
law [Eq. iJSJl] 



a = a a cxpHTm/T) 1 / 4 ], 



(57) 



with parameter Tm increasing with I according to Tw oc 
l/.9s oc (Z// s ) /3+2 . The growing Tm leads to exponen- 
tially increasing resistivity, represented by the ascending 
branch of the curve in Fig. [3 Physically, the suppression 
of the DC conductivity stems from decreasing density 
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of low-energy states available for transport, just like in 
conventional doped semiconductors^ or commensurate 
CDW systems^ 



VI. CONCLUSIONS AND COMPARISON WITH 
EXPERIMENT 

It is widely recognized that interactions must play a 
significant role in determining the properties of ID and 
quasi- ID conductors, because in such materials the di- 
mensionless strength of the Coulomb interaction is very 
large, r s 3> 1. In the presence of impurities, these sys- 
tems behave as insulators and do not possess metallic 
screening. Thus, the interactions are both strong and 
long-range. Our main goal in this paper was to under- 
stand the effect of such interactions on the nature of the 
low-energy charge excitations and their Ohmic dc trans- 
port. To that end we formulated a generic model of an 
anisotropic electron system with strong Coulomb interac- 
tions and disorder and presented its theoretical analysis. 
We elucidated the origin of the low-energy charge excita- 
tions in this model and demonstrated that their density of 
states possesses a soft Coulomb gap. In 3D case, we found 
that the Coulomb gap exhibits a power-law dependence 
on the energy distance from the Fermi level. We dis- 
cussed how the prefactor and the exponents of this power- 
law vary as a function of the impurity concentration and 
other parameters of the model. We also discussed how 
the Coulomb gap is manifested in the variable-range hop- 
ping conductivity at low temperatures. 

One of the central results of our theory is a nonmono- 
tonic dependence of a on the impurity concentration N, 
as shown in Fig. |2 where we sketched a as a function 
of I = 1/Na\, i.e., as ./V decreases, at fixed T. As clear 
from that Figure and the discussion in Sec. [HJ at large 
N the conductivity increases as N decreases, similar to 
behavior found in metals. In contrast, at small N the 
conductivity drops as N goes down, which resembles the 
behavior of doped semiconductors. For intermediate N, 
our theory predicts the existence of an A-independent 
plateau. 

Another way to represent these theoretical predictions, 
common in semiconductor physics, is shown in Fig. El In 
that Figure the dependence of the logarithm of conduc- 
tivity on the inverse temperature is depicted for a series 
of samples, each with fixed N. An unusual circumstance 
illustrated by Fig. is the crossing of the curves that 
correspond to different samples. In Fig. [HJl we show that 
up to two crossing points may exist between one curve 
that corresponds to I < l s (curve 2) and another curve for 
I > l s (curve 4). The higher-T crossing occurs when the 
curve 4 goes through the activation regime, the lower- 
T one — when it exhibits the Mott VRH. The dirtier 
(/ < l s ) sample obeys the ES behavior in both instances. 
Another property we tried to emphasize in Fig. is the 
role of the 2/5-law as the upper bound of the conduc- 
tivity regardless of the the sample purity. For samples 
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FIG. 6: Logarithm of the conductivity vs. the inverse temper- 
ature for several samples labelled in the order of increasing I, 
i.e., sample purity, (a) Curve 1 corresponds to I — h < l s and 
displays only the ES regime. The higher-T activation regime 
is beyond the limits of the graph. Curve 2 is for I = l 2 > l 3 , 
and so the both the activation behavior and the ES law are 
visible. Curve 3 is for I = £3 > I2. It shows the complete 
sequence of the transport regimes: the activation, Mott, 2/5, 
and ES laws. Curve 4 depicts the behavior of a sample with 
an impurity concentration another notch lower than that of 
sample 3. In panel (a) it goes through the activation and the 
Mott regimes. In panel (b) that covers considerably lower T, 
curve 4 also exhibits the 2/5 law, followed by the ES law. The 
2/5-law is shown by the solid line along the lower edge of the 
shaded region. The curves can skim along this line but cannot 
cross it. The higher the sample purity the lower the tempera- 
ture at which the sample starts to exhibit the 2/5-law but also 
the wider the range of T over which this law persists. Curve 
5 depicts the case / 3> h where any kind of VRH transport 
corresponds to Ohmic resistances higher than the experimen- 
tal measurement limit, so that only the activated transport 
can be observed. 



with low impurity concentration the 2/5-law is also the 
envelope curve, see Fig. [SJd. 

Let us now turn to the experimental situation. The 
transport behavior of a number of organic compounds, 
including TMTSF-DMTCNQ, TTF-TCNQ, and NMe- 
4-MePy (TCNQ) 2 is indeed in a qualitative agreement 
with our theory. It should be clarified that such materi- 
als form CDW phases that in addition to the usual 2k p 
periodicity, have an appreciable or, in some cases, even 
predominant Akp harmonics. The latter is considered the 
evidence for the strong Coulomb interaction^ and so is 
precisely the case we studied in this paper. 

In the experiments of Zuppiroli et alm^L the trans- 
port in the aforementioned compounds was studied as 
a function of defect concentration, which was varied 
in situ by irradiation of samples by high-energy parti- 
cles. Admittedly, the nature of the such defects is not 
known with certainty. Some suggestions in the litera- 
ture include atomic displacements, broken bonds, poly- 
mer cross-linking, and charged radicals. At the same 
time, the effect of the irradiation on transport seem not 
to depend much on the type of particles used (X-rays, 
neutrons, or electrons) and instead to correlate primarily 
with the total absorbed energy^ This fact is interpreted 
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TABLE II: Fitting parameters for data displayed in Fig. [7| 



Quoted defect concentration, c 


Tes, K 


E A , K 


-6% 




900 


3.5% 




380 


2.3% 


2800 




1.9% 


3700 





as evidence that microscopically defferent defects influ- 
ence the transport in electron crystals in a similar way, 
so long as they act as strong localized pinning centers. 
Under this assumption, it is legitimate to compare the 
data from the irradiation experiments with our theory 
even though so far we assumed that defects are created 
by charged acceptors (see Sec.Pl. We do so in some detail 
below. 

In Fig. \7\ we show an extensive set of data on trans- 
port in irradiated TMTSF-DMTCNQ that we assembled 
by digitizing Fig. 2 of a review article by Zuppiroli— and 
original references therein. Apparently, some data series 
in this figure represent the same sample with successively 
increasing radiation dose, and some correspond to phys- 
ically different specimens. For simplicity, we refer to all 
of them as different samples. The percentage labels on 
the plot are the estimates of the molar concentrations of 
defects (c) quoted by the experimentalists. The points 
on the c = trace are from an unirradiated sample. 

As shown in Fig. [7] the data for the two most disor- 
dered samples can be successfully fitted to the activation 
law and the next two samples — to the ES formula. This 
transition from the activation to the ES law with increas- 
ing disorder is in agreement with our theory (cf. curves 
1 and 2 in Fig. EJ) . The obtained fit parameters Ea and 
Tes are given in Tabled Both Ea and Tes scale roughly 
linear with c, in agreement with Eqs. (|llf) and (O- From 
Eq. IjlQI) we deduce that kIc — na ~ 1 nm, which has 
the correct order of magnitude (assuming k *~ 1). One 
should keep in mind here that the absolute numbers for 
c were obtained by the authors of Ref. OH] using certain 
arguable assumptions. In our opinion, the scaling with c 
may be more reliable than the absolute values quoted be- 
cause (if no annealing occurs) the relative magnitude of c 
should scale linearly with the irradiation time, known to 
experimentalists without any fitting parameters. Com- 
bining Eqs. ijjjj and Ijllfl we further deduce that a±/£± 
and l/^x ratios are some modest numbers less than ten, 
as may be expected from Eqs. <|25|) and (JSJl. 

As a final remark on high-disorder samples, we would 
like to mention that the scaling of the longitudinal dilec- 
tric constant k x with the defect concentration (irradi- 
ation time) consistent with Eq. 1221) was reported in a 
separate set of experiments on Qn(TCNQ) by Janossy 
et alMi Together with the transport data, this makes 
a compelling case for the validity of the metallic-rod 
(interrupted-strand) model for organic electron crystals 
with I <C l s - For such systems we can claim a semi- 
quantitative agreement with the experiment. 




0.1 0.2 

i/r(K-') 



FIG. 7: Low-temperature conductivity of TMTSF-DMTCNQ 
samples damaged by X-ray radiation. The dots on c > 1.9% 
and c < 0.35% curves were generated by digitizing the ex- 
perimental data in Fig. 1 of Ref. |^ and Fig. 2 of Ref. I3BI 
respectively. The percentages stand for the defect concentra- 
tions quoted in those papers. The solid line through the 6% 
and 3.5% data are the best fits to the activation; through the 
2.3% and 1.9% — to the ES law. In both cases the prefactors 
(a a and 00) were taken to be T- independent. The thin line 
through the 0% data (unirradiated sample) is the best fit to 
the Mott law based on the T < 20 K points. 



Let us now turn to the conductivity of weakly damaged 
samples. 35 As one can see from Fig. they show metal- 
lic behavior at high-T, a conductivity maximum at the 
Peierls temperature of about 42 K, and a decrease in er, 
i.e., the semiconducting behavior, at lower T. As T drops 
by a factor of two or so below the Peierls temperature, 
the decrease of a with 1/T becomes considerably more 
gentle than an activation law. One can not help notic- 
ing a similarity between the behavior of c = 1.9,0.35, 
and 0% samples in Fig. [7| and curves 1, 3, and 4, re- 
spectively, in Fig. There is also an unambiguous ev- 
idence for the existence of the crossing point between, 
e.g., c = 0% and c « 1.2% traces at T = 21 K (see be- 
low). However, an attempt to fit the c = 0% data to the 
Mott law is not particularly successful, see Fig.0 There- 
fore, we only wish to emphasize a qualitative agreement 
with our prediction that for a fixed T, the conductiv- 
ity of "clean" and "dirty" samples should show oppo- 
site trends, see Figs. |21 and EJ Indeed, the conductivity 
of the low-disorder samples (c < 0.35%) increases with 
the radiation dose in contrast to the behavior shown by 
the high-disorder samples (c > 1.9%) where it decreases. 
In fact, another experiment showed this contrasting be- 
havior in a great detail. In that experiment^ In a was 
measured at the fixed temperature of T = 21 K, while 
defect concentration was varied essentially continuously 
over the range of < c < 2.5%. It was found that a 
initially increases by two orders of magnitude, reaches 
a maximum, and then drops by five orders of magni- 
tude as a function of c. Overall, this is in a qualitative 
agreement with Fig. |2 except instead of the well-defined 
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2/5-plateau, In a shows only a broad maximum. Similar 
features are demonstrated also by TTF-TCNQ and NMe- 
4-MePy (TCNQ) 2 , see Fig. 1 in Ref.ES, and Figs. 11 and 
12 in Ref. [3^. We conclude that for low-disorder CDW 
organics our theory agrees with experiment in some gross 
qualitative features. The quantitative agreement can- 
not be verified because the dynamical range of measured 
conductivities is too narrow. Further low-temperature 
experiments are desired to clarify the situation and to 
prove or disprove the existence of the novel 2/5-law. 
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Let us now switch to inorganic CDW. Several com- 
ments are in order. The electron-electron interactions 
in these materials are also very strong, 10 ^ r s ~ 100. 
However, inorganic CDW are predominantly 2kp, and 
there is an ample evidence for the important role of the 
electron-phonon coupling in the CDW dynamics. This 
coupling can lead to an enhancement of the electron ef- 
fective massjifi^ that would result in a short localiza- 
tion length. If the mass enhancement is indeed large, the 
VRH transport should be observable only in materials 
with short hopping distances, i.e., large impurity con- 
centrations. Examples include highly doped bronzes^* 
and perhaps, the Pt-chain compound KCP&S From 
this perspective, the reports of a VRH-like transport in 
a relatively pure samples of TaS3 [Refs. and|3| and 
blue bronze [Refs. and 0] are puzzling and require 
further investigations. 



Finally, let us comment on another broad class of 
quasi-lD systems, conducting polymers. In comparison 
to CDW, polymers have a much higher degree of struc- 
tural disorder and a complex morphology that depends 
on the preparation method. Typical samples contain of 
a mixture of crystalline and amorphous regions, with the 
correlation length^ of the order of 10 nm. In the un- 
doped state polymers are commensurate CDW semicon- 
ductors with a Peierls-Mott energy gap 2 " ~ 1 eV. Doping 
shifts these systems away from commensurability point 
and suppresses the gap but it is often inhomogeneous and 
is an additional source of disorder. At low and moder- 
ate doping the T-dependence of the conductivity often 
resembles ES and/or Mott VRH laws, see a short review 
in Ref. 44. A systematic study of the VRH conductivity 
dependence on doping has been attempted by Aleshin et 
al~ In those experiments doping of PEDT/PSS samples 
was varied by controlling the pH of the solution at the 
sample preparation stage. It was observed that at pH < 4 
the VRH exponent A [Eq. (|T7|) ] was close to 0.5 and TVrh 
decreased with pH, while at larger pH, A was close to 0.4 
= 2/5 and Tvrh did not depend on pH. This resembles 
the behavior that follows from our theory, provided the 
concentration of the pinning centers decreases with pH. 
We leave the tasks of extending our theory to the case 
of conducting polymers and explaning these intriguing 
experimental results for future investigations. 



APPENDIX A: DIELECTRIC FUNCTION AND 
ELASTICITY IN A CLEAN QUASI-1D 
ELECTRON CRYSTAL 

In this section we derive expressions for the elastic 
moduli and the dielectric function of a pure crystal. 

Following the literature on interacting ID electron sys- 
tems^ and CDW. 48 we describe dynamics of electrons 
on ith chain by a bosonic phase field tfi(x,t). The the 
long-wavelength components of electron density n, is re- 
lated to if i by rii = d x (fi/2TT. The elasticity theory of 
the system can be formulated by identifying the elastic 
displacement u with (2ir/a)(p and taking the continuum 
limit. Neglecting weak interchain tunneling and dynam- 
ical effects we choose our starting effective Hamiltonian 
H in the form 

H = fdx[H +H c ], (Al) 
Ho = \C° X "I + E J H co <^ ~ W)> ( A2 ) 

i ij 

He = dxdx 1 rii(x)Uij (x — x')rij(x'), (A3) 

ij 

Let us briefly describe the notations here. The Hamilto- 
nian is split into the short-range (H ) and the long-range 
Coulomb (Hq) parts. Jy represents the Coulomb cou- 
pling between the CDW modulations of electron density 
on chains i and j. C° is the charge compressibility of a 
single isolated chain. In a large-r s ID system C° is dom- 
inated by the exchange-correlation effects, leading tc 62 

2e 2 a 

C£~-— In | (A4) 
k R 

where R is the characteristic radius of the electron charge 
form- factor in the transverse direction, i.e., the "radius" 
of the chain. In most physical realizations, we expect 
R <C a and negative C®- The positivity of the elastic 
modulus of the system, required for thermodynamic sta- 
bility, is recovered once we take into account the long- 
range part He of the interaction energy, parametrized 
by the kernels f/y. f/y is defined as the bare Coulomb 
kernel Uo(r) = e 2 /nr convoluted with the single-chain 
form-factors F, e.g., -F(q_i_) = exp(— q\R 2 ). 

We are interested in a linear response where the cosines 
in Eq. i|A2|l can be expanded in ip thereupon the effec- 
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tive Hamiltonian becomes quadratic. If an external elec- 
trostatic potential V oxt (q) acts on the system, the total 
equilibrium potential V(q) will in general contain Fourier 
harmonics with wavevectors q + G where G are the re- 
ciprocal vectors of the 2D lattice formed by the trans- 
verse coordinates of the centers of the chains. We define 
the dielectric function e(q) of the system as the ratio 
t^,xt(q)/t^(q) for q in the Brillouin zone of this lattice. 
Via standard algebraic manipulations in the reciprocal 
space we arrive at 



e(q) 



S*(q) = C° x 



Aire q. 



2 _2 



1 



4 q 2 B x (q)ql+B±qj 



E 



F 2 (G) 



Ka l <^ ll + (<U + G ) 2 ' 



Bj_ = 4 7 r 2 a 2 L ^J 4: 



(A5) 



(A6) 



(A7) 



In the limit q x <C a 1 , q± -C a ± we obtain Eq. l|3t)|) 
reproduced below for convenience, 



e(q) = k - 



q 2 q 2 + aq]_ ' 



(A8) 



Here ro = {C x k/ Aire^^ax is the Thomas-Fermi 
screening radius, C x > is the effective charge compress- 
ibility of the system, 



C x — C x + 



47re 2 v-^ F 2 (G) 2e 2 n x 

Ka, K a 



and a = B±/C x is the dimensionless anisotropy param- 
eter. Let us now discuss some consequences of Eq. (|A8ll . 

(i) The dielectric function has the same form as in 
quasi-lD systems with small r s (see, e.g., Ref . E?^ . 

(ii) The screening radius rr> is of the order of the 
interchain separation ax- Therefore, Brazovskii and 
Matveeenko's results 2 ^ for the soliton energy W and the 
soliton length l s remain qualitatively correct for high 
r s provided we use td ~ a ±- 

(iii) The interaction energy U(r) of two point-like test 
charges separated by a large distance r can be calculated 
by the Fourier inversion of 



Aire 2 
e(q)9 2 



(A10) 



For ax 



U(r) 



> r 2 D one finds (cf . Ref. |2 



■ exp 



1 



2«|a;| V 2r D ^/ax 2 + r]_ + y/a\x\ J 

(AH) 

The potential U (r) is not exponentially small only within 
the paraboloid r 2 ^ < y/ar D \x\ [cf. Eq. (J20J)]. At the sur- 
face of such a paraboloid we have r± -C ^/a|a;| and U(r) 
acquires a simpler form quoted in Sec. 



U(r) 



2k\x 



exp 



iy/a re | 



(A12) 



(iv) Finally, the effective longitudinal and transverse 
elastic moduli of the system are given by 



Y x = 



C x 
a 2 a\ 



9 2 ' 



Y ± =aY x . 



(A13) 
(A14) 



APPENDIX B: TUNNELING IN A 3D CRYSTAL 
WITH A LOW IMPURITY CONCENTRATION 

The localization lengths £j_ and £ x needed for calcu- 
lation of the VRH transport are determined by long- 
distance tunneling of charge excitations. The problem 
of tunneling is nontrivial because a broad spectrum of 
charge excitations exists. Leaving more detailed inves- 
tigations for future work, we concentrate on two possi- 
ble tunneling mechanisms: by electron- like quasiparticles 
and by many-body excitations, the 2-7r-solitions. 

In the quasiparticle mechanism the charge is carried 
by a single electron while all other electrons remain un- 
perturbed in their quantum ground states. The ratio- 
nal for examining this mechanism is its minimal possible 
tunneling mass. The problem of calculating £ x and £j_ 
reduces to the quantum mechanics of a single particle in 
a fixed external potential. Clearly, the optimal tunneling 
path should go through the interstitial positions where 
the energy barrier is the lowest. It is convenient then to 
formulate the problem as a problem on a lattice of such 
interstitial positions. The relevant variables are the on- 
site energies and the hopping matrix elements. The on- 
site energies are all equal to e; n t ~ e 2 / Ka. The hopping 
terms for the interchain tunneling, tx, are determined by 
the band-structure in the case of tunneling inside a chem- 
ically synthesized materials. In the case of tunneling be- 
tween distant ID conductors t± cx exp(— ax/as)- The 
hopping matrix clement for tunneling along the chain, 
t x , can be estimated straightforwardly, with the result 
t x = £int exp(— CV-y/rJ), CV ~ 1. In the absence of impu- 
rities, the problem is reduced to the propagation through 
a periodic lattice. The eigenstates in that model are la- 
belled by wavevectors k, according to the tight-binding 
dispersion relation 

Et6(k) = £ int -2t a: cos(A; a: a_L)-2tx[cos(fc !/ ax)+cos(fc z a_L)]. 

Below the band edge e = £j nt — 2t x — 4tj_, the eigen- 
states are exponentially decaying with distance. The 
corresponding localization (decay) lengths can be related 
to the imaginary parts of the complex k solution of the 
equation £tb(k) = e. In the case of interest, e -C £int; 
tx , t x <C £i„t , we obtain 



a = 
?± = 



\n(e int /2t x ) 
ax 

]n(e int /2txY 



a 



(Bl) 
(B2) 



where the superscript q stands for "quasiparticle." Al- 
though Eq. (|B2|) was derived for a clean system, it is clear 
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that impurity do not affect this result unless present in 
gigantic concentrations (I ~ a). Indeed, an individual 
impurity can modify the local on-site energy by at most 
a numerical factor, while £x(x) depend on the on-site en- 
ergy only weakly, logarithmically. In principle, impu- 
rity clusters with atypically low on-site energies, resonant 
with e, do exist but as well known from the analysis of 
the resonant tunneling problem in random systems, such 
events are exponentially rare and do not contribute to 
the bulk localization length in any appreciable manner. 

Let us now turn to the solition mechanism, we we at- 
tempt to profit from the fact that in the bulk, the soliton 
is the charge excitation of the lowest possible energy, so 
that the energy barrier could perhaps be lower and the 
tunneling more effective. However, in the case of inter- 
chain tunneling, this is not the case. Indeed, the direct 
tunneling of a soliton to a different chain is impossible 
because the soliton is a composite many-body excitation. 
The closest to the interchain soliton tunneling that one 
can imagine is a two-stage process, where, one electron 
first tunnels to the adjacent chain, and then, on the sec- 
ond stage, it pushes away other electrons in the region of 
length l s to form a soliton. Since the initial energy bar- 
rier is still eint and the charge spreading only increases 
the tunneling action, it is clear that such a contrived 
process offers no advantage compared to the simple one- 
stage quasiparticle mechanism. Therefore, £j_ is deter- 
mined by the latter and coincides with leading to 
Eq. (JSJ. Note that for the case of distant chains where 
t± ~ exp(— oi/ob), the correct limiting result £j_ = as 
is recovered. 

In contrast, for the tunneling along the chain the soli- 
ton mechanism is the winner. Consider first the longitu- 
dinal tunneling of a soliton in the absence of impurities. 
Empoying the usual imaginary-time picture, the action 
for tunneling over a distance x ^> l s can be estimated as 
the product of the energy barrier W and the tunneling 
time r ~ x/u. Here it, given by 



ma 



1/2 



1 



(B3) 



is the sound velocity*** The tunneling amplitude is of the 
order of exp(— Wx/hu). Using W ~ e 2 /nl s , we arrive at 
the estimate of the localization length as follows 

I, 



(i = 00), 



(B4) 



where the superscript s stands for "soliton." Clearly, 
S ^ £x' so that the soliton mechanism dominates the 
longitudinal tunneling. To account for the dilute impuri- 
ties, we should add to the above expression for the action 
an extra term (x/l)[hy/r^ln(l s /a)}. Here the factor (x/l) 
is the average number of impurities on the tunneling path 
of length x and the expression inside the square brack- 
ets is the action cost for compactification of the charge-e 
from the length l s to length a and spreading it back dur- 
ing the tunneling through each impurity. In this manner, 



we obtain a corrected expression for £ x , which coincides 
with Eq. (J5SJ. 



APPENDIX C: DIMENSIONAL ENERGY 
ESTIMATES FOR THE COLLECTIVE PINNING 

In this Appendix we use the ideas of collective pinning 
to derive the growth of elastic distortions in a quasi-lD 
crystal pinned by strong dilute impurities. We also derive 
the estimates of the corresponding gain in the pinning 
energy density. 

We start by reformulating the argument leading to 
Eq. (|33fl in the language conventional in the literature 
devoted to collective pinning^ To do so we note that 
since the energy of a given soliton dipole E s depends on 
the background elastic displacement field u, each impu- 
rity exerts a force / = —dE s /du ~ W/a on the crystal. 
The long-range variations of u appear in response to such 
random forces. Let Ait(D) be a characteristic variation 
of u over a distance D in the transverse direction and let 
X be a typical distance over which a variation of the same 
order in the x-direction builds up. Our next step is to 
estimate the total energy E of a volume V — X x D x D 
(relative to the energy of a pristine crystal) . 

The energy consists of elastic, Coulomb, and pinning 
parts, 



E = E cl + E c + E p 



(CI) 



In its turn, E el is the sum of the longitudinal and trans- 
verse terms, 



E e 



V. 



(C2) 



The Coulomb energy is of the order of U(q x ,q±)p 2 V, 
where U(q) — 47re 2 / 'nq 2 is the Coulomb kernel, q x — 
1/X, q± = 1/D are the characteristic wavevectors in- 
volved, p = end x u ~ enAu/X is the charge density asso- 
ciated with the londitudinal compression, and n = \ j aa\ 
is the average electron concentration. Below we show 
that D <C X, so that q±_ » q x , U(q) ~ e 2 D 2 /n, and 
finally, 



E u ~ — — 
k X 



Au 

nn 2 



(C3) 



The pinning energy can be estimated as E pin ~ 
—AuJ2j fj: where fj ~ W/a is the force exerted on the 
lattice by jth impurity. The average number of impuri- 
ties in the volume V is Ni = NV and fj have random 
signs; hence, E pm ~ — (W/a)Au v / W l , or 



E p " 



-W 



Au 
a 



D 

ax 



1/2 



(C4) 
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Combining Eqs. (|A13|I and HC2IC4|) . we arrive at 



W f Au 

a \a± 

'Au\ ( D 



aX 



D 4 
Xa\ 



a 



1/2 



(C5) 



X and Au can now be found by optimizing E for a fixed 
D. Not surprisingly, we find that X and D are related by 
the defining equation (|30[1 of the paraboloid introduced 
in Sec. IV Al (for rn ~ a± case), 



X 



D 1 



(C6) 



Such a paraboloid is an invariable feature of the elas- 
tic response of the quasi-lD crystal to external forces. 
What is surprising however is that Am 2 ~ {l s /l)a 2 is 
small and does not depend on D, at odds with Eq. (|32[) . 
The resolution of this contradiction comes from a real- 
ization that what Au(D) really represents is the elas- 
tic distortion due to the adjustment of the crystal on a 
single scale D. In fact, there is a hierarchy of smaller 



scales D, D/2, D/A, . . . , r™ ln , on which adjustments are 
approximately independent. The correct estimate of Au, 
Eq. (|32[1 . is obtained once we sum over all such scales, 
Au 2 ~ Ma 2 l s /l, where M ~ ln(D/r^ in ) is the number 
of scales. Thereby, we recover Eq. (|32ll and as an addi- 
tional benefit, we find the expression for the energy of 
the collective pinning, 



e ■ -ii i — y (— 2 

a J Vax. 



(C7) 



Let us define the pinning energy density by £p- m — 
E/XD 2 . Using Eqs. PH |) . (35 | . £3), and JUT)) we ob- 
tain the estimate of £ pm at the Larkin scale as given by 
Eq. I]34p. As clear from this derivation, both Eq. and 
(|34|l are essentially the lowest-order perturbation theory 
results. It is generally expected that the growth of Au 2 
with r slows down beyond the Larkin length 65 and that 
adjustments on larger scales do not lead to any substan- 
tial increase in the pinning energy density. In this case, 
Eq. (|34fl is the final estimate of £ p [ n in the thermody- 
namics limit. 
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